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1. Introduction 

The atmospheric neutrino results reported by the Super-Kamiokande [1] 
and MACRO [2, 3] experiments indicate that neutrinos oscillate, these data 
being consistent with z/,j «-» v T oscillations. The small value of the difference 
of the squared masses (5 x l(T 4 eV 2 < Am 2 < 6 x KT 3 eV 2 ) and the 
strong mixing angle (sin 2 29 > 0.82) suggest that these neutrinos are nearly 
equal in mass as predicted by many models of particle physics beyond 
the standard model. Also, the LSND experiment [4] supports z/„ <-> v e 
oscillations (Am 2 < 0.2eV 2 ); some solar neutrino experiments [5] suggest 
that u e could oscillate to a sterile neutrino v e <-> v s (Am 2 ~ 10 _5 eV 2 ). 
The direct implication of neutrino oscillations is the existence of a non-zero 
neutrino mass in the eV range or lower, and consequently a not negligible 
hot dark matter (HDM) contribution ^ to the total energy density 
of the universe. Recent works [6, 7, 8] show that the addition of a certain 
fraction of HDM component to the total energy density of the universe can 
lead to the agreement between the Cosmic Microwave Background (CMB) 
anisotropy power spectrum at the small scales and the observations of the 
Large Scale Structure (LSS) of the universe as derived by the redshift galaxy 
surveys [6]. 

The CMB anisotropy pattern contains information related to the phys- 
ical processes occurring before the last scattering of the CMB photons; 
the LSS data reflects the clustering regime effects in our local universe. 
These measurements can provide independent probes for the structure of 
the universe on similar comoving scales at different cosmological epochs. 
This paper discusses the implications of massive neutrinos properties for 
CMB and LSS measurements. 
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2. Properties of the neutrino background 

2.1. COSMOLOGICAL LIMITS ON NON-DEGENERATE NEUTRINO 
MASSES 

Evidences has been accumulated that our universe is a low matter den- 
sity universe (see e.g [9] and the references therein). Indications like the 
Hubble diagram of Type la supernovae and the acoustic peak distribu- 
tion in the CMB anisotropy power spectra (see e.g. [10] and the references 
therein) point to a universe dominated by vacuum energy (cosmological 
constant A) that keeps the universe close to flat. The combined analysis 
of the latest CMB anisotropy data and Type la supernovae data indi- 
cates fi m = 0.25lg];| and Q\ = 0-63^Q 23 [10] f° r matter and vacuum 
energy densities respectively, inferring a Hubble constant value of Ho ~ 
65 km s- 1 Mpc" 1 (h = H /W0 km s^Mpc" 1 = 0.65). Adding a HDM 
component to the ACDM model (ACHDM) leads to a worse fit to LSS 
data, resulting in a limit on the total non-degenerated neutrino mass of 
m v < 2eV for a primordial scale invariant power spectrum and m v < 4eV 
for a primordial scale free power spectrum [8]. A stronger upper limit is 
obtained from the matching condition of the LSS power spectrum normal- 
ization eg (defined as the rms amplitude of the galaxy power spectrum in 
a sphere of radius 8/i -1 Mpc) at the COBE scale and at the cluster scale 
[9]. For the case of the ACHDM model having £l m = 0.3, U A = 0.7 and 
a primordial scale invariant power spectrum it is found an upper limit of 
m v < 0.6 if Hq < 80 km s -1 Mpc -1 . Also, the constraints on the cosmolog- 
ical parameters obtained by using the most recent CMB anisotropy data 
[11] combined with Type la supernovae data implies a best fit model with 
£l m ss 0.33, 1^=0.67 and a neutrino density parameter Q u « 0.12 (J2i m Vi 
« 92h 2 n u = 4.66eV) if priors H =65 km s" 1 Mpc" 1 and ft 6 h 2 =0.02 are 
assumed [12]. 

2.2. COSMOLOGICAL LIMITS ON NEUTRINO DEGENERACY 

A way to improve the agreement of the cosmological models involving neu- 
trinos with the higher values of Hubble parameter is to consider that neu- 
trinos are degenerated. The neutrino degeneracy parameter is defined 
as £„ = fi u /T u , where fi u is the neutrino chemical potential and T v is the 
neutrino temperature. The neutrino degeneracy enhances the contribution 
of the HDM component to the total energy density of the universe [13] 
changing the neutrino decoupling temperature [14], the abundance of light 
elements at the big bang nucleosynthesis (BBN) [15], the CMB anisotropies 
and the matter power spectrum [16]. Constraints on neutrino degeneracy 
coming from BBN [15] indicate —0.06 < £„ e < 1.1 and |^ | < 6.9. Bounds 
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on the neutrino degeneracy parameter are obtained by combining the cur- 
rent observations of the CMB anisotropies and LSS data [16]. The analysis 
of the CMB anisotropy data obtained by the Boomerang experiment indi- 
cates bounds on massless neutrino degeneracy parameter of jT | < 3.7, 
if only one massless neutrino species is degenerated, and \£ u \ < 2.4, if the 
asymmetry is equally shared among three massless species [17]. 



2.3. MASSIVE NEUTRINO MIXING AND THE PHASE SPACE 
DISTRIBUTION FUNCTION 

In the standard assumption of the mixing of massive neutrinos, the neutrino 
flavor eigenstates are described by a superposition of the mass eigenstate 
components that propagate differently with different energies, momenta and 
masses [21]. This paper considers two massive neutrino flavors and v T 
and a third neutrino flavor u e is massless. The mixing occurring in vacuum 
between and v T can be written as: 



v lt ) = cos 6>o | V2 )+ sin O \u ?! ) (1) 
= - sin^ol^} + cos^ol^), 



\v. 



where V2 and v% are the mass eigenstate components and 6q is the vacuum 
mixing angle. It is usual to consider that each neutrino of a definite flavor 
is dominantly one mass eigenstate [21]. In this circumstance, the dominant 
mass eigenstate component of is V2, that of v T is vj, and their difference 
of the squared masses is Am 2 = — m 2 .. 

In the expanding universe, the system of neutrinos drops out from ther- 
mal equilibrium with e ± , photons, the small fraction of baryons and other 
massive species (e.g. pairs) when the ratio of the averaged weak in- 

teraction rate to the expansion rate falls below unity. The present values 
of neutrino temperature T UQ and the present temperature of the CMB pho- 
tons T ss 2.725K are related through T UQ /T ~ (3.9/g*/)) 1 / 3 , where g* D 
is the number of degrees of freedom in equilibrium at neutrino decoupling. 
For a degeneracy parameter < 15 the neutrino decoupling temperature 
is of few MeV [14] and from the entropy conservation one obtains that 
after decoupling T VQ ~ (4/ll) 1 / 3 To. It follows that the total degeneracy pa- 
rameter £ v is conserved and the present lepton asymmetry of the neutrino 
background is [15]: 

T — _ ' -^o tc3 l J^t 1 
"~ 12C(3) T [ ^ + *" J ' 

where £(3) is the Riemann function of 3. At decoupling, neutrino with 
mass in eV range behave like a relativistic particle and its full phase space 
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distribution function can be written as a pure Fermi-Dirac distribution 
function f®(q), a zero-order distribution that depends only on momenta, 
plus a perturbed part \I/(x, q, a) that depends on momenta, position and 
time [22]: 

/„(x, q, a) = + *(x, q, a)) f v {q) = - 1 , (2) 

where: E v = \Jq 2 + a' 2 m'f, is the neutrino energy, a is the cosmological 
scale factor (ao = 1 today), q is the comoving neutrino momentum, ±£„ is 
the neutrino(-sign)/antineutrino(+sign) degeneracy parameter and vn v is 
the dominant mass component. The full neutrino phase space distribution 
function was computed through numerical simulations [23] based on the 
standard particle-mesh (PM) method. The initial neutrino positions and 
velocities was generated from the HDM matter density fluctuation power 
spectrum by using Zel'dovich approximation [24]. Figure 1 shows (left pan- 
els) the contour plots of the constant particle probabilities in Sq-SpnDM 
plane, where Sq is the fluctuation of the neutrino comoving momentum and 
Sphdm is the fluctuation of the neutrino density field. It is also shown (right 
panels) the momentum distribution functions obtained from numerical sim- 
ulations (continuous line) compared with a pure Fermi-Dirac distribution 
function (dashed line) . The time evolution of the phase space distributions 
f v and f Ur of 1*^ and v T can be written as [25]: 

8f u ^,q,t)/dt = -<P>, M /, M (x,q,t)+<P>, T /, T (x,q,t), 
df UT (x,q,t)/dt = -<V> Vt /, T (x,q,t)+<P>^/^(x,q,t), (3) 

where: < V > v and < V > Ut are the neutrino mixing probabilities, dt = 
da/(aH) and H is the Hubble expansion rate: 

o 

In the above equation G is the gravitational constant, f2 m = Qf, + Q c dm + 
is the matter energy density parameter and Qcdm: are the energy 
density parameters for baryons, cold dark matter and neutrinos, Q r is the 
energy density parameter for radiation (including photons and one rela- 
tivistic neutrino flavor), is the vacuum energy density parameter and 
flk = 1— £l m — is the energy density parameter related to the curvature of 
the universe. The set of equations (3) assume that, apart from the neutrino 
mixing, no other physical processes take place, leading to the conservation 
of total neutrino number densities n u and degeneracy parameters £ u . It 
follows that non- vanishing mixing and degeneracy parameters would lead 
to changes of the energy density, pressure, their perturbations, the shear 
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Figure 1. Left panels: the contour plots of constant particle probabilities in the 
5q — 5pHDM plane. From exterior to interior the contours correspond to: 0.75, 0.5 and 0.25 
probability. Right panels: the full neutrino phase space distribution functions obtained 
from numerical simulations (continuous line) and the pure Fermi-Dirac distribution func- 
tion (dashed line). 

stress and the energy flux of each massive neutrino flavor, quantities that 
are related to the neutrino phase space distribution function [25]. Figure 
2 presents the time evolution of the energy density parameters and of the 
density perturbations of different components (matter, radiation, vacuum, 
neutrinos) for few values of neutrino mixing and degeneracy parameters. 
One can observe that the variation of Am 2 and A£j, affect the energy den- 
sity parameters and the density perturbations for all components, while the 
variation of sin 2 29q affect only the neutrino density perturbations, reflect- 
ing the conservation conditions imposed by equations (3). 

3. CMB temperature and matter density fluctuations 

During the epochs where the anisotropics are generated, neutrinos with 
masses in eV range can have significant interactions with the photons, 
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Figure 2. The evolution with the scale factor of the energy density parameters f2; 
and of the density perturbations Si of the various components in a ACHDM model hav- 
ing Qt = 0.023, Qa = 0.67, Q m = 0.33, h = 0.65, ra v = m Vfl + m VT — 0.6eV and 
£» = + £„ T = 5. (a) and (al): sin 2 26 = 0.8, A£„ = 4, Am 2 = 0.24eV 2 (continu- 
ous line) and Am 2 = 0.18eV 2 (dashed line), (b) and (bl): Am 2 = 0.24eV 2 , A£„ = 4, 
sin 2 2^0 = 1 (continuous line) and sin 2 29o = 0.45 (dashed line). Panels (c) and (cl): 
Am 2 = 0.24eV 2 , sin 2 28 = 0.45, A£„ = 4 (continuous line) and A£„ = 2 (dashed line). 



baryons and cold dark matter particles only via gravity. Thus differences 
introduced by the neutrino mixing of the phase space distributions and the 
lepton asymmetry of the neutrino background are due to the differences 
in the gravitational field of the neutrinos. As it was shown in the previous 
section, the neutrino energy density and pressure as well as the perturbed 
the energy density and pressure, the shear stress and energy flux are modi- 
fied by the mixing of the neutrino degenerated background. They determine 
differences in the equation of state of massive neutrinos (w = p u /Pu), in the 
expansion law of the universe (eq. 4) and originate the metric perturbations 
affecting the growth of fluctuations in other components. 
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Figure 3. (a): The dependence of the time evolution of neutrino free streaming distance 
and sound horizon distance on Am 2 and A£„ when: Am 2 = 0.24eV 2 and A£„ = 4 
= 0.16] (continuous line), Am 2 = 0.2eV 2 and A£„ = 4 [Q„ = 0.15] (dashed line), 
Am 2 = 0.24eV 2 and A£„ = 1 [fi„ = 0.07] (dash-dotted line). The vertical line indicates 
the approximated value of the scale factor when massive neutrinos start to make the 
transition from radiation to matter, (b): The growth functions of the perturbations at the 
present time (z — 0) obtained for: Am 2 = 0.24eV 2 and A£„ = 4 [fi„ = 0.16] (continuous 
line), Am 2 = 0.2eV 2 and = 4 [Q„ = 0.15] (dashed line), Am 2 = 0.24eV 2 and 

A£„ = 1 = 0.07] (dash-dotted line). The cosmological model is a flat ACHDM model 
having: Q b = 0.023, fi A = 0.67, Q, m = 0.33, h = 0.65, m v = m„„ + m Vr = 0.6eV and 

£,v = + £,v T = 5. 

3.1. THE NEUTRINO GRAVITATIONAL INFALL 

Characteristic to the cosmological models involving massive neutrinos is 
the scale-dependence of the growth rate of perturbations. Neutrinos can 
cluster via gravitational instabilities only on distances below a character- 
istic distance, the free streaming distance, analogous to the Jeans scale of 
self-graviting systems (see e.g. [22] and the references therein). 
The neutrino free streaming distance is related to the causal horizon dis- 
tance 77(a) through [26]: 



1 77(a) 




da 



(5) 



kf s y/1 + (a/a nr ) 



a 2 H 



where: kf s is the free streaming wave number, a is the cosmological scale 
factor, H is the Hubble expansion rate given by the equation (4) and a nr 
is the scale factor when massive neutrinos start to become non-relativistic 
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Figure 4- Maxima of the mass distribution of the cold dark matter component (panel 
a), baryonic component (panel b) and hot dark matter component (panel c) obtained 
from numerical simulations of 32 3 particles in a box of size 32 Mpc for the ACHDM 
model [fl b = 0.023, fl A = 0.67, Q m = 0.33, h = 0.65, m v = m„„ + m Vr = 0.6eV and 
£,v = £ Vft + £,v T = 5]. Panel d): The mass distribution of the HDM component obtained 
from numerical simulations for different neutrino fractions (see also the text). The total 
mass in all simulations is M = 2.27 x 1O 15 M . 



(a nr = (1 + z nr ) 1 ss 3k B T U0 /m u c 2 ). 

Another important scale is given by the sound horizon distance of the 
baryon-photon fluid defined as [26]: 

rv(a) „ 11 
R s = c s (a)dr}, c s (a) = -- (6) 

Jo 6 1 + I a HT 

where: c s (a) is the sound speed and f2 7 and are the energy density pa- 
rameters of photons and baryons at the present time. Panel a) in Figure 
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3 presents the dependence of time evolution of the neutrino free streaming 
distance and sound horizon distance on Am 2 and A£„. Panel b) of the same 
figure shows the evolution with the wave number of the growth functions 
of perturbations obtained for the same values of these parameters. The in- 
crease of Am 2 when A£ u is kept constant leads to the increase of Q. v and 
of z nr . The same effect is obtained when Atti 2 is kept constant and A^ is 
increased. At early times, when neutrinos are relativistic, the free streaming 
distance is approximately the sound horizon distance. After the neutrinos 
become non-relativistic, the free streaming distance decreases with time 
becoming smaller than the sound horizon distance. Also, the growth of 
perturbations is suppressed at smaller scales, the magnitude of the sup- 
pression depending on neutrino parameters. 

It follows that neutrinos can cluster gravitationally on increasingly small 
scales at latter times, damping the amplitude of the density perturbations. 
This implies that the characteristics of the cosmological structures at small 
redshift, as the mass density profile and the velocity dispersion profile of 
clusters, can probe the neutrino properties. Figure 4 presents maxima of 
the mass distribution of the cold dark matter component, the baryonic 
component, and the hot dark matter component, obtained from numeri- 
cal simulations in a flat ACHDM model. The cold dark matter and the 
baryonic matter components accrete the hod dark matter component via 
gravity. The gravitational potential governing this process is given by the 
Poisson equation V 2< £ = 4TrGa 2 pe m (k,a), where p = £l m Pcr is the averaged 
density, p cr is the critical density and e m (k, a) is the net perturbation of the 
density field. Also, panel d) of the Figure 4 presents the mass distribution 
of the HDM component (obtained for each simulation bin and normalized 
to the total mass of each bin) for different neutrino fractions. 

3.2. CMB ANISOTROPY AND THE MATTER POWER SPECTRA 

In the linear perturbation theory, the CMB anisotropy and matter transfer 
function are computed by the integration of coupled and linearized Ein- 
stein, Boltzmann and fluid equations [22] that describe the time evolution 
of the metric perturbations in the perturbed density field and the time evo- 
lution of the density fields for all the relevant particles. Figure 5 presents 
the dependence of the CMB anisotropy power spectra (left panels) and 
matter density fluctuations power spectra (right panels) on the neutrino 
difference of the squared masses Am 2 , the mixing angle sin 2 26q and the 
lepton asymmetry L v [25]. For all cases it is assumed a total neutrino 
mass m u = m v + m UT = 0.6eV and a total neutrino degeneracy parameter 
= £ Vfi +£ VT = 5. The CMB power spectra are normalized to COBE/DMR 
four- Year data and the matter density fluctuations power spectra are nor- 
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Figure 5. The CMB anisotropy power spectra and the matter density fluctuations power 
spectra. Panels a) and al): sin 2 28o = 0.8, A£„ = 4 and (from top to bottom) Am 2 =0.24, 
0.2, 0.01 eV 2 . Panels b) and bl):Am 2 = 0.24cV 2 , A£„ = 4 and (from top to bottom) 
sin 2 2#o = l, 0.8, 0.45. Panels c) and cl):Am 2 = 0.24eV 2 , sin 2 26> =0.8 and (from top 
to bottom) A£„=4, 3, 1. In each panel are presented (dashed lines) the power spectra 
obtained in the same cosmological model when Am 2 = 0, A£„ = and sin 2 2Go = 0. 



malized to the abundance of the rich clusters at the present time [23] . 
The CMB anisotropy power spectra show two distinct features: a verti- 
cal shift of the C\ at large I (present in all panels) with the increasing of 
Am 2 , sin 2 26*o and A£„ values, that results from the modification of neu- 
trino free streaming scales, and an horizontal shift of the Doppler peaks 
to lower I when Am 2 and A£ u are increased, that results from the in- 
crease of the sound horizon at recombination. The horizontal shift of C\ 
is not present in panel bl) because the variation of sin 2 2#o, when Am 2 
and A£ u are fixed, does not change Q u and consequently the Hubble ex- 
pansion rate is not changed. The matter density fluctuations power spec- 
tra reflects the same features: the increasing of Am 2 and A^ u values re- 
sults in the modification of neutrino free streaming scales and consequently 
a suppression of the growth of fluctuations on all scales below the free 
streaming scale (kf s ~ 0.02Mpc _1 ). By using the Fisher information ma- 
trix method it was found [25] that the future high precision CMB experi- 
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ments as Planck surveyor will permit to detect a lepton asymmetry of the 
neutrino background h u > 4.5 x 10 -3 , a difference of the neutrino squared 
masses Am 2 > 9.7 x 10 _3 eV 2 and a mixing angle sin 2 29q > 0.13. Also, the 
combined data from Planck and Sloan Digital Sky Survey will permit the 
determination of the total neutrino mass m v > 3 x 10 _2 eV, with any prior 
information on cosmological parameters [23]. 

4. Conclusions 

The standard Cold Dark Matter model has difficulties in matching the CMB 
and LSS measurements. Cosmological models involving a mixture of CDM 
and HDM particles, the CHDM models, are able to fit, the excess large 
scale power seen in galaxy surveys and the CMB temperature fluctuations. 

The results from neutrino oscillation experiments indicate a non-zero 
neutrino mass in the eV range, or lower. This implies a non-negligible con- 
tribution of neutrinos to the total energy density of the universe. The ex- 
perimental evidence indicating a present low matter density universe, dom- 
inated by the vacuum energy (Q m ~ 0.3, f^A — 0.7) and a higher Hubble 
parameter value (Hq ~ 65 km s _1 Mpc _1 ) are in agreement with cosmolog- 
ical models involving neutrinos if one considers the lepton asymmetry of 
the neutrino background. 

At the times when the anisotropics were generated, neutrinos had signif- 
icant interactions with the photons, baryons and cold dark matter particles 
only via gravity. Neutrinos can not cluster via gravitational effects on scales 
k below the free streaming scale k < kf s (large wavelengths). For k > kf s 
the growth of the density perturbations is suppressed, the magnitude of 
this suppression depending on Q m , Q u , Am 2 and A£„. 

The neutrino homogeneous quantities (density and pressure) and inho- 
mogeneous quantities (density and pressure perturbations, the shear stress 
and the energy flux) are also changed in the presence of non-degenerated 
neutrino mixing, leaving imprints on the CMB temperature fluctuations 
and the matter density fluctuations power spectra. 

The determination of the fundamental cosmological parameters with 
high precision CMB and LSS surveys requires complementary knowledge 
of neutrino properties from cosmic rays and long base-line experiments. 

Acknowledgements: I acknowledge the organizers and V. Berezinsky 
for useful discussions during this workshop. 
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